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ABSTRACT: 

In this paper are presented the results of numerical computations of 

the effect of space charge upon the image size and image distance in a 

converging beam of charged particles. The minimum beam radius (image 

size) and its position (image distance) are plotted as a function of beam 

“13 “6 

perveance ranging from 10 to 10 . The beam convergence angle was 

varied from about 0.03 to 0. 15 radians and calculations were carried out 

for 3 values of aperture radius (0.1 cm, 0.25 cm and 0.50 cm). The 

, . r , 1 TT + 7 + 20 + , 8 5 . + 

results are given for electrons, H , Li , Ne , and Rb . 
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I. INTRODUCTION 

In recent years there has been an increasing demand for low energy 
ion beams, well defined in space and energy, to be used in experimental 
investigations of atomic collision processes. In the energy region below 
100 eV the effects of space charge upon usable beam intensity, beam shape 
and ion energy become so serious as to represent the limiting factors in 
many experiments . 

Pierce has designed a widely used lens utilizing curved electrodes 
to create the fields required to simulate flow between concentric spheres . 
Although the Pierce gun has been used for the production of low energy 
ion beams with some success, it was recently shown by Simpson and 

O 

Kuyatt that there exists a fundamental limitation on the beam intensity 
which can be obtained by any unipotential gun operating at low energies. 
They show that in order to attain the current densities required for 
saturation, in the presence of thermal effects at the emitter surface, 
one would have to exceed the space charge limitations imposed by the 
diode nature of any unipotential gun. From this it is concluded that the 
range of applicability of unipotential guns lies at comparatively high 
energies and wide beam convergence angles. To operate outside these 
limits it is necessary to decouple the thermal energy and space charge 
effects by first accelerating the particles to a fairly high energy and 
then decelerating them to their final low energy. Simpson and Kuyatt^ 
have designed an electron gun using this multistaging principle and we 
have adapted their design for use with low energy ion beams. 

A schematic diagram of the gun is given in Fig. 1 . The first stage 
extracts the particles at high energies and forms a beam which is then 
decelerated to the desired final energy in the second stage. The first, 
or extraction stage is based on designs given by Soa^ and the decel- 
eration stage is based on the properties of plane electrodes, as given 
by Spangenberg and Field. ^ The electrode potentials V and V a are 
measured with respect to the emitter which is maintained at the desired 
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FIGURE 1 
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beam potential, V q , while the exit electrodes are grounded. 

In addition to circumventing the basic problem described above, this 

design has considerable mechanical advantages over the conventional 

Pierce gun. All electrodes are plane; displacement of the deceleration 

stage with respect to the extraction stage and variation of the electrode 

potentials allows one to obtain a wide variety of beam shapes, intensities 

and energies quickly and conveniently. 

The details and results of our experimental investigation of this gun 
0 

are reported elsewhere and will only be summarized very briefly here. To 
obtain a low energy ion beam we used a thermal emitter (L^O'A^O^ *n SiC^) 

as a source of Li ions . Current densities of 10 to 10 amps/cm^ at 

beam energies of 2 to 100 eV were obtained and both the spatial and energy 

distributions of the ions were investigated. Minimum beam radii of 1 to 2 

mm occurred a few centimeters beyond the last element of the gun. Beam 

profiles were measured as a function of axial distance and well collimated 

beams of up to 16cm were obtained. The energy spread in the beam was 

found to be about 0.22 eV in good agreement with the expected thermal 

energy spread from the source which operates at about 1200°C. Mass 

7 + 

analysis shows that beam purities are well in excess of 99% Li when 
isotopically enriched emitter material is used. 

II. BEAM DYNAMICS 

The beam geometry in the absence of space charge effects is shown 

in Fig. 2. As the beam leaves the aperture of radius R the rays ema- 

B 

nating from the edge are initially directed as shown and those from 

interior points are contained within the outer rays so as to be focused 

toward a spot of radius r^ a distance L away. The origin of coordinates 

is such that the beam travels along the z axis, the coordinates of the 

aperture is z and of the image z . The beam is assumed to be axially 
d a 

symmetric. 

In considering the space charge effects we shall follow the work 
of Glaser and the analyses of Watson and Wendt . We shall assume 
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FIGURE 2 
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a beam of charged particles of sufficient density to be considered a 
continuous medium for purposes of constructing energy functions. 

The effect of the space charge of the beam is to deflect the rays 
away from the axis. We assume that the current density is uniform 
across the beam cross section and zero outside the outer rays, and that 
the total current is constant. If the convergence of the rays is suffi- 
ciently small then a small segment of the beam forms a cylinder of 
nearly constant radius in which the axial forces are nearly zero, and 
hence the axial velocity z is very nearly constant and given by 

A _ / 2-e.V ' 

^ V w\ 



( 1 ) 



where e = charge on the particle, V = beam potential, and m = particle 
mass . 



We shall now proceed to calculate the trajectories of the outer 
and inner rays as shown in Fig. 2 under the influence of space charge 
(in field free space) subject to the simplifying assumptions given above. 



Trajectory of the Outer Ray 

The radial component of the electric field due to the space charge 
can be written 

i 

l 



YA 



2-eV 

where i is the total current, eV is the beam energy, and r^ is the radius 



r, 

it, eV is 

of the outer ray at any point. Let 

- -Sr 

be the charge to mass ratio of the particles. The radial equation of 
motion of a particle on the outer ray is then 

/ e. t 



( 2 ) 



wr, * 



If we define the constant K. to be 

1 « 

K, = " 
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we may write Eq. (3) as V 2- 

2k_ 

i n 



(3) 
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(5) 



Multiplying both sides by dr^ and integrating gives 

(£) 

where the constant of integration r is chosen such that at 
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r, * n 



r = (O 



(7) 
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Using equations (1) and (4) we can write j 

n--i (%)]'■ 

where we choose the negative root in order to restrict consideration to the 
portion of the beam between the aperture and the f ocal point, i.e. , the 

region where the beam is convergent and r < 0. This gives on integration 

^ - 1 









T’j 



where the primes denote variables of integration. 



Finally, setting 



= 1a (^)l 






we find 



where 



= - 2K 



E 
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(8) 



(9) 






If we define the Dawson function 

D&> = j e** 

then Eq. (9) becomes 

= - 2 K [l>GO - £ G$e>] 
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( 10 ) 



(11) 



Since tables of D(S) are available we can evaluate equation 11 in 
terms of the constant K. This equation gives r^ as an implicit function 



of z and the parameters R r , r^,L, T^and the beam perveance P. 

To exp 
perveance *P 



B* a' 

To express K in terms of the minimum radius r o and the beam 

l 



we note from above , that 



K = K 



V. 

K, 



( 12 ) 
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where 



Y - 

* ' V z\y 



Thus 



K - r c / mJB- 



or 






K = Cv 0 t 
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l 

where 



c - 



(4) 



(13) 



Trajectory of the Inner Ray 

We know from Gauss' Law that only the charge enclosed by the 

inner ray will affect its trajectory. Assuming a uniform radial charge 

distribution, the ratio of the current included within r 9 to the current 

2 Z 

included within r-, will be r 2 and since the field depends only on the 
enclosed charge we can use Eqs . (2) and (4) to write 



^ - K, -±- 

and the radial equation of motion of a particle on the inner ray is 

•• Y" 

where r^ is given by 

r i - 4 ex ? ^**1 

from Eq . (8) . 

We may differentiate parametrically as 

5 gW} _ ^ ^ S 






els d r, dS 



r, 



where f is given by (6). ds/dr, may be found from (8), giving 
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( I ^ M *lr. 



and 
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(14) 



(15) 
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Putting (18) in (15) gives 

_ *) * d r*. 



__ . - 25 - 2r u * O 

_ ds 



( 19 ) 



which has the solution-*- ^ 

W = C, {s 1 ^ [> +■ C t 9<C s ->3 

where and C 2 are constants of integration, and 

m W" £ 644 ^ 0,4 " 

0 (s) is the well-known error function which is extensively tabulated'*''*- . 

We evaluate the constants of integration by noting that at z = z^ 



( 20 ) 



(21) 



(see Fig . 2) 



and 
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( 22 ) 



Note also that from equation (9) 

d? 

Applying Leibnitz' rule gives 



tion (9) 1 

377 ' - 2 K ) e**> js' 
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and 



whence 



From (16) 
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dr, 
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(4tt) b = ~ 
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r 0 u 
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so 



>i • ■Es e*e {- 



Now from (20) 

exip ^5^ L 



1-hC 






and 



Hi" 



— 



= C, 



? l s *i at ir& 5 ^ +C2 -^ +c ^ J 
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(25) 

(26) 
(27) 



Applying Leibnitd'rule to (21) gives 

- -S— exp ^-s*^ 



(28) 



and 



<As VtF” 

a* - 



i 



2.S K 



(29) 



Putting (23), (28), (29), and the second equation of (22) into (27) and 
evaluating at ^ = Rg / gives 

where 0g denotes 0 (Sg). Solving (26) for and putting this and (24) in 
(30) gives r 

2 fr k:^ B r K 



(30) 



i + c i. 4>% 



r.» L 



Let 



u = a ^ 

With this we have immediately 

lvl 



^ - 



and from (26) ^ ^ rQ 

C-, -V-, (> - M <z( B ) 

so the complete solution for ^ (see equation 20) is given by 

K = K> e *p jj" 



(31) 

(32) 

(33) 

(34) 
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where r is given by (25), N by (31), s by (24) and <b by (21). It should 

o B 

be noted that equation (34) gives r as a function of the constants R , 

^ B 

r , L, r\ and the beam perveance P at every value of z. 

0. 



Calculation of Minimum Beam Radius and Postion 

The spreading of the rays under the influence of space charge is 

illustrated in Fig. 3. We shall distinguish between two cases depending 

on whether the value of the beam perveance P is larger than or smaller 

than a critical perveance P defined below: 

cr 

Case 1; P>P (Fig. 3a): For this case the outer ray reaches 
cr 

its minimum before its intersection with the inner ray 

originating at the diametrically opposite side of the 

aperture. Thus the outer ray reaches its minimum 

value r, = r before -r. = r.. and hence the minimum 
1 0 2 1 

radius occurs at z = z . 



Let z -z_ = L_ . At r =r , s = 0 (from equation 8). 

° B F 1 ° 

Since D(o) = 0 we get from equation (11) 

Lp ' D(sg^ ( 35 ) 

and r o = r^ is given by equation (25). 

Case 2; P < P cr (Fig* 3)5): Here the outer ray will be intersected 

by the inner ray originating at the diametrically opposite 

side of the aperture, before r^ reaches its minimum value 

r . If we call this point of intersection z„ then at 
0 r 

Z = Z F <• “ r 2 = r i = V 



101 Z F ' Z B = h • 



r l " r 2 



The minimum radius occurs where 
Equation (16) and (34) gives 



i- N (y B - <fi r ) =-i 



from which 






Z. 

Nl 



(3 6) 
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( 37 ) 



F 

and 



can be found from (36) and we find Tp from 
= r Q exp | s p 2 | 

= 2 K [D (s b ) - P (s F )] 



(38) 



When P = P cr , Lp (case 1) = Lp (case 2) and the method 
of either case may be employed in evaluating Tp and Lp . 



III. RESULTS OF NUMERICAL COMPUTATIONS 
In this section we present the results of numerical computations for 
the image distance Lp (Equ. 35 and 38) and the image radius Tp (Equ. 25 
and 37) . Each of the figures (Fig. 4-18)^ shows 3 curves of Lp as a 
function of beam perveance "P ( ^ for three values of 
the beam convergence angle y (Fig „ 2) (tan y - 0 . 03 , 0.09 and 0.15). 

On the same figure we plot Tp as a function of P for the same three 
values of tan y . The curves are aligned so that values of Lp and rp 
can be read for any given value of P and the same value ©f tan y by 
drawing a vertical line at the proper value of perveance. Three values 
of the final aperture radius Rp (See Fig. 2) (Rp = 0.1, 0.25 and 0.50 
cm) and five mass values were used: electrons, ^-H + , ^Li + , ^®Ne + 
and 85 Rb + . 
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